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The temperature dependence of Coulomb blockade peak 
height correlation is used to investigate how adding electrons 
to a quantum dot alters or "scrambles" its electronic spec- 
trum. Deviations from finite-temperature random matrix the- 
ory with an unchanging spectrum indicate spectral scrambling 
after a small number of electrons are added. Enhanced peak- 
to-peak correlations at low temperature are observed. Peak 
height statistics show similar behavior in several dot configu- 
rations despite significant differences in correlations. 
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Electron transport through irregular quantum dots — 
i.e. micron-scale islands of confined charge weakly con- 
nected to electronic reservoirs - are expected to, and in 
some cases actually do, exhibit universal statistics associ- 
ated with quantum chaos JjJ. An example where theory 
and experiment agree well is the distribution of Coulomb 
blockade (CB) peak heights [0-p|. At temperatures T 
that are much smaller than the mean level spacing of the 
dot A, transport on a CB peak is mediated by resonant 
tunneling through a single level HQ . Large fluctuations 
in CB peak heights in this regime reflect the fluctuat- 
ing strength of coupling of the chaotic wavefunction in 
the dot to the modes in the leads, leading to universal 
statistics sensitive only to time-reversal symmetry [^],[|, 
in good agreement with experiment ^,^]. 

At higher temperatures, A < kriT < Ec, where Ec is 
the classical charging energy, each CB peak contains con- 
tributions from ~ fc^T/A quantum levels, and one would 
expect roughly this number of adjacent peaks to be cor- 
related in height. This assumes that the spectrum of the 
dot does not change as electrons are added. On the other 
hand, if adding electrons alters the spectrum, then the 
correlation length in peak number, n c , will not grow be- 
yond a certain value, m, which roughly measures (but is 
not equivalent to) the number of added electrons needed 
to completely "scramble" the electronic spectrum. 

This Letter presents measurements of the temperature 
dependence of the CB peak-to-peak height correlation 
and peak height statistics for gate-confined GaAs quan- 
tum dots, and compares these results to finite tempera- 
ture random matrix theory (RMT) calculations that ne- 



glect spectral scrambling ||. We find that the number 
of correlated peaks n c (T) saturates at m ~ 2 — 5, with 
smaller dots saturating at smaller m. At the low tem- 
perature end, we find that n c (T) is larger than the value 
predicted by RMT. That is, correlations in peak heights 
beyond thermal smearing exist for reasons that are not 
clear. Some possible explanations are considered below. 

In contrast to the dependence on dot configuration 
found for the peak height correlations (as reflected in 
n c (T) and m), peak height statistics are found to be very 
similar for all device configurations. This suggests that 
peak statistics are more robustly "universal" than peak 
correlations, not surprising considering that distributions 
are not sensitive to spectral scrambling. The ratio of the 
standard deviation to mean of peak heights is found to 
be smaller than predicted, possibly due to the effects of 
decoherence. 

What does one expect to be the effect of adding elec- 
trons on the spectrum of a quantum dot? In the limit of 
weak electron-electron interactions (and neglecting shape 
deformations caused by changing gate voltages) a fixed 
spectrum of single-particle states is simply filled one at 
a time, leading to n c (T) ~ fc^T/A and m ^> 1. In the 
opposite limit of strong interactions, the spectrum could 
be totally scrambled with the addition of each electron, 
giving m ~ 1. For a GaAs quantum dot containing many 
(~ 100 or more) electrons, RPA calculations j|,[l0| (ap- 
propriate for weak interactions) indicate that fluctuations 
in the ground state energy due to interactions are small, 
of order r s g _1 / 2 A where g is the dimensionless conduc- 
tance of the dot and r s is the so-called gas parameter, 
the ratio of potential to kinetic energy of the electrons 
(r s ~ 1 — 2 in GaAs heterostructures). For a ballistic- 
chaotic dot containing M electrons g ~ A/" 1 / 2 , giving a 
rough estimate for the number of electrons needed to 
scramble the spectrum, 1 < m <^ A/" 1 / 2 /?" 2 , assuming 
that fluctuations accumulate randomly as electrons are 
added to the dot. 

Measurements of CB peak spacing statistics have in 
some cases found rms fluctuations in Ec as large as 15% 
[0,[L2|, consistent with classical estimates |13 and nu- 



merics |llj for strong interactions (where RPA fails), 
suggesting that one or a few electrons can significantly 
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alter the charge arrangement of the dot. Other experi- 
ments |0| have found smaller fluctuations, of order A, 
suggesting a lesser degree of rearrangement. In all ex- 
periments the peak-spacing distribution is found to be 
roughly Gaussian, which is surprising considering that 
spin degeneracy would naively imply a bimodal distribu- 
tion. Suggested explanations include wave-function de- 
pendent interactions |^,[l5| and scrambling due to shape 
deformation p6| . Magnetofingerprints of CB peaks in 
dots with ~ 50 — 100 electrons show a persistence of 
both the addition and excitation spectrum over at least 
6 peaks, with changes in the excitation spectrum (mea- 
sured by nonlinear magnetotransport) of order A upon 
adding one electron, leading occasionally to an exchange 
of a pair of levels in the spectrum 0. Symmetric, few- 
electron dots also show a sort of scrambling in the sense 
that the ground state shell-filling structure may differ 
from the corresponding excited state before the N th elec- 
tron is added, for as few as TV = 4 electrons [Q. For 
such small systems (TV <~ 5 — 6) exact calculations are 
possible, and a statistical approach to scrambling is not 
necessary. 

Before describing the experiment, we discuss the gen- 
eralization of the theory of CB peak height fluctuations 
to temperatures that are comparable or greater than A, 
but without including spectral scrambling Well- 
formed CB peaks require weak tunneling from each of 
the levels A coupling left and right leads to the dot, with 
tunneling rates (r^,T^) <C A/ft,, and also low bias and 
temperature, (eVd s , ksT) <C Ec, where Vd s is the voltage 
bias across the dot. In this regime, the peak conductance 
G r , 



has the form 
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where a(T) = ^ A a\W\(T) is a weighted sum of normal- 
ized lead-dot-lead conductances, a A = 2rfr^/(r(r, A + 
T^)). The dot is assumed to be symmetrically cou- 
pled to the leads, with average tunneling rates Tf — 
r, A = r/2. In the experimentally relevant regime 
(fceT, A) <C Ec, the thermal weights w\(T) are given by 
w x (T) = 4/(AF^ - E F ){n x ) M [l - f(E x - E F )], where 
AFaa = Faa — Fjv"-i is the difference in the canonical free 
energy of TV and TV" — 1 non-interacting electrons on the 
dot, {n\)j^ is the canonical occupation of level A with TV 
electrons on the dot, Ep = [Ep + er]V g — (Af—l/2)Ec] is 
the effective Fermi energy with V g tuned between TV — 1 
and TV" electrons on the dot, E\ is the energy of level A, 
and /(e) = 1/(1 + e e / keT ) is the Fermi function. Equa- 
tion (^) generalizes previous results for low temperatures, 
hT <C fee? 1 ^ A, and yields the known distributions for 
a in that limit 

Within a noninteracting model, no correlations be- 
tween neighboring peak heights are expected for ksT <C 
A. At higher temperatures, correlations appear as each 



level is able to contribute to several nearby peaks. For 
both numerical RMT data and experimental data, we 
compute a discrete correlation function C(n) from a se- 
quence of TV peaks (G max )i 



C(n) = 
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where 8g i = {{G max )i - (G max ) NjB ) is the fluctuation 
of the i th peak height around the average (calculated 
over both peak number and magnetic field). The cor- 
relation length n c is then calculated from a Gaussian fit, 
C(n) — e~(™/" c )" . The gaussian form is not based on any 
theoretical model but appears to accurately describe the 
shape of both the RMT and experimental data. 

The quantum dots we discuss are fabricated using e- 
beam lithography to pattern Cr/Au gates on the sur- 
face of a GaAs/AlGaAs heterostructure 900 A above the 
2DEG layer. Multi-gate dot design allows changes in 
device size and shape by changing gate voltages. Ac- 
tive control of point contact gates during sweeps of the 
"plunger" gate V g compensates any unintentional capaci- 
tive coupling, allowing many peaks to be swept over with- 
out changing the average transmission of the leads. All 
data were taken using two-wire lock-in techniques with 
2 iieV bias at 13 Hz. An electron base temperature of 
45 mK and ratio r\ of gate voltage to dot energy were 
extracted from the CB peak width versus temperature 
using standard methods . The charging energy Ec in 
each configuration was calculated using r] and the average 
peak spacing in gate voltage JtJ . The mean level spacing 
A was measured from the sub-structure (corresponding 
to excited states) in the differential conductance at finite 
bias using a small ac signal added to a dc bias. Measure- 
ments were made with broken time-reversal symmetry, 
with 3 — 10 </> through the device area. For each temper- 
ature, V g was scanned over 50-100 peaks and B was then 
changed by ~ <fi /Ad t to give independent peak height 
statistics. 

Figure |l| shows two typical series of CB peaks at lower 
and higher T, where the increased correlation between 
peaks at higher T is clearly evident. Peak height fluctu- 
ations 6g i (right insets in Fig. |]) are extracted from each 
series using fits to cosh" 2 lineshapes around each peak, 
and the correlations C(n) are then calculated according 
to Eq. (||). Plots of C(n) for two dots are shown in Figs. 
||(a, b). For both configurations, the correlation length n c 
increases with increasing T at lower T . However, for the 
smaller configuration n c saturates above 300 mK, while 
for the larger configuration n c continues to increase well 
above this temperature. In still larger devices, saturation 
is not observed up to temperatures of 700 mK. 

Correlation lengths n c (T) for three measured dot con- 
figurations (A = 20, 28, and 38 ^eV) are shown in Fig. 
||(c) along with the RMT results. Each data point in 
Fig. E|(c) represents data from ~ 500 CB peaks. The 
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FIG. 1. (a,b) Coulomb blockade peaks in conductance g as 
a function of gate voltage V g at (a) 45 mK and (b) 400 mK 
from device 1. Insets: SEM micrograph of device 1. Peak 
height fluctuations Sgi extracted from these data sets. 



RMT curve was computed by applying Eq. to a peak 
sequence that is similar in length to the experimental 
data, and generated according to Eq. ([!]) assuming a 
uniformly-spaced spectrum |8) . The RMT results do not 
change significantly when Wigner-Dyson statistics for the 
spectrum is included. The saturation of n c (T) at m ~ 2 
for ksT > 0.5 A for the smallest device is evident in Fig. 
H(c). The larger dot begins to saturate for larger n, with 
a larger ratio ksT/A, suggesting that the spectrum of 
the larger dot is less prone to scrambling. The observed 
scale of saturation, m, as well as the trend for m to in- 
crease with Af, appears consistent with the RPA estimate 
given above. 

As the gate voltage is swept, two distinct changes oc- 
cur, both of which can cause spectral scrambling. The 
first is that the number of electrons and size of the dot 
change; the second is that the shape of the dot changes 
due to local movement of the boundary at the position 
of the gate. This second effect was considered recently 
in Ref. [ [l6| to explain the nearly-Gaussian peak spacing 
distribution seen in several experiments |Tl],[l2|,[l4j] . The 
two effects can be separated using a dot with two plunger 
gates, which allows pure shape distortion without chang- 
ing M by increasing one gate voltage and decreasing the 
other. In practice, it is easier to raster over the two gate 
voltages, as seen in Fig. ||(a). Horizontal and vertical 
directions correspond to single-gate CB measurements, 
while a downward diagonal following a single peak cor- 
responds to pure shape distortion with fixed J\f. Corre- 
lations in the same dot measured at fixed N (measured 
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FIG. 2. (a,b) Peak height correlations C(n) at 45 mK, 
100 mK, 200 mK, 300 mK and 400 mK for (a) dot 1 and 
(b) dot 2. (c) Temperature dependence of correlation length 
n c for different device configurations, and numerical RMT 
result. Inset: Grayscale plots of conductance for 3 successive 
CB peaks, showing paired peaks i and i + 1, presumably a 
spin pair. 



along diagonals) and changing M (measured along hor- 
izontals) can be compared by evaluating both correla- 
tions in terms of V g \ rather than n. Comparing C(AV g i) 
for the two cases (Fig. ||(b)) shows that the correlation 
length associated with shape deformation is larger by a 
factor of ~ 4 than that associated with a changing N '. 
This indicates that the saturation of n c (scrambling) is 
dominated by changes in electron number rather than by 
shape distortion. Further work is needed to determine if 
this result is universally true, but it appears to hold in a 
variety of gate-confined dots that we have measured. 

All dot configurations show an enhanced correlation 
length n c (T) at low T compared to RMT, as seen in 
Fig. ||(c). Thus the data suggest that temperature alone 
does not explain the enhancement, at least within RMT. 
Several possible explanations for this have been proposed 
]l5| , ^9|p0| ; we further note that correlations at low T can 
result from similar peak heights in spin-paired levels (see 
Fig. ||(c), inset) which may appear in either adjacent or 
non-adjacent peaks, depending on the size of interaction- 
induced spin splitting. 

Finally, we investigate the statistics of peak height fluc- 
tuations as a function of temperature. It is convenient 
to consider a normalized distribution to remove any tem- 
perature dependence of the average peak height, (a(T)). 
We define a — a/ (a), so that the distribution P{a) 
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FIG. 3. (a) Grayscale conductance plot of CB peaks as 
a function of V g i and V 3 2 from device shown in inset of (b), 
with A ~ 23 fieV, at 90 mK. The appearance of peaks as 
short vertical bars reflects the coarser sampling of V g 2 com- 
pared to V g i ; the patterns of connected bars moving up and to 
the right are an artifact of this display, (b) Correlation func- 
tion C(AV" 9 i) of peak height fluctuations for fixed N (dashed 
curve) and fixed V g 2 (changing Af, solid curve). 



will have (a(T)) = 1 for all T. The standard deviation 
<j{a) = y/ (a 2 ) — 1, which characterizes the width of the 
distribution relative to its mean (i.e., cr(a) = a(a)/ {a}), 
is shown for three dot configurations along with RMT re- 
sults in Fig. ^. All experimental data show very similar 
temperature dependences despite significant differences 
in correlations described above. We conclude that peak 
height distributions, which are not sensitive to scram- 
bling, show more "universal" behavior than correlations. 
Notice, however, that the experimental data all have 
smaller height fluctuations than predicted by RMT. This 
can also be seen in the full distributions P{a) comparing 
experiment (histograms) and RMT (solid lines), shown as 
insets in Fig. ^|. The departure from finite-temperature 
RMT is likely due to decoherence effects, and might pro- 
vide a novel tool for measuring decoherence in nearly 
isolated structures. Ongoing work in this direction is in 
progress. 
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FIG. 4. Temperature dependence of scaled peak height 
distribution width, a (a), along with numerical RMT result 
(solid curve). Insets: Full distributions of scaled peak heights 
P(a) for all dot configurations combined, at feT = 0.1 A, 
k B T = 0.5 A, and k B T = 1.5 A, along with RMT distribu- 
tions. RMT results use a uniformly-spaced spectrum (rather 
than Wigner-Dyson) . 
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